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Abstract: We use the GALPROP code and the Advanced Composition Explorer (ACE) data to derive 
the cosmic ray (CR) isotopic composition at the sources. The composition is derived for two propagation 
models, diffusive reacceleration and plain diffusion. We show that the compositions derived assuming 
these two propagation models are different. We also compare the isotopic composition at the sources 
with the latest solar composition. 
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Introduction 

CR source abundances are normally derived using 
the leaky-box or weighted slab approximations to 
interstellar propagation. While this is a valid pro- 
cedure under some conditions (see [20] for a re- 
view) at least for stable nuclei, there are reasons 
for preferring a more physically-based approach. 
For example, the distribution of CR sources (e.g., 
supernova remnants) in the Galaxy is probably 
peaked towards the inner Galaxy, as recently sup- 
ported by 7-ray observations [19], and this af- 
fects the path-length distribution in a way that re- 
quires a detailed spatial propagation model. An- 
other example is the effect of diffusive reaccelera- 
tion, which is probably important in reproducing 
the energy-dependence of the secondary/primary 
ratios, and this will also affect the derivation of 
source abundances. The parameters of the models 
used to derive the source abundances should also 
be compatible with other observational constraints 
from gas surveys, 7-rays, synchrotron and so on 
[17, 12, 18, 19]. 

In this paper, we use the GALPROP CR propaga- 
tion code [16] to derive source isotopic abundances 
from ACE data. Both the astrophysical model and 
the cross-sections play a key role in the uncertain- 
ties in such a computation; here we present pre- 



liminary results, reserving detailed discussion to a 
future journal paper. 

IMethod 

The GALPROP code' computes a complete net- 
work of primary, secondary and tertiary produc- 
tion starting from input source abundances, as de- 
scribed in [16, 18, 14]. The nuclear reaction 
network is built using the Nuclear Data Sheets. 
The isotopic cross-section database is built us- 
ing the extensive T16 Los Alamos compilation of 
the cross-sections [8] and modern nuclear codes 
CEM2k and LAQGSM [9]. The most important 
isotopic production cross-sections (2H, 3H, 3He, 
Li, Be, B, AI, CI, Sc, Ti, V, and Mn) are calcu- 
lated using our fits to major production channels 
[11, 12]. Other cross-sections are calculated us- 
ing phenomenological approximations by Webber 
et al. [21] (code WNEWTR.FOR versions of 1993 
and 2003) and/or Silberberg and Tsao [15] (code 
YIELDX.Ol lOOO.FOR version of 2000) renormal- 
ized to the data where it exists. The K-capture and 
electron stripping processes are included, where a 
nucleus with one electron is considered a separate 
species because of the difference in lifetime. 



1. http://galprop.stanford.edu 
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Figure 1 : Quality of the fit: fractional deviations of propagated elemental abundances from ACE observa- 
tions [22]. 



The propagation equation is solved numerically 
starting with the heaviest nucleus (i.e., ^*Ni), com- 
puting all the resulting secondary source functions, 
and then processes the nuclei with A—1. The pro- 
cedure is repeated down to A = 1. To account for 
some special /3^-decay cases (e.g., ^°Be— ^^"B) the 
whole loop is repeated twice. The current version 
employs a full three-dimensional spatial grid for all 
CR species, but for the purposes of this study the 
2D cylindrically symmetrical option is sufficient. 

To calculate the isotopic source abundance, we 
adopt an iterative procedure, which uses the devi- 
ations of calculated propagated abundances from 
abundances observed by ACE at 200 MeV/nucleon 
[22] to correct the source abundances until a good 
fit is obtained. In practice very good agreement 
(a few %) can be obtained in about ten iterations. 
The basic parameters of the propagation model 
are based on fitting the energy-dependence of the 
B/C ratio, for which the best data are available 
and cross-sections best known. The parameters are 
given in [14]. Solar modulation is calculated using 
the force-field approximation [6]. The modulation 
potential $ = 450 MV corresponds approximately 
to the period of solar activity when the data were 
collected. An analysis of source abundances of iso- 
topes of P, S, Ar, Ca is given in [13]. Discussion 
of propagation of Li, Be, B in the same models can 
be found in [4] . 



Results 

We applied the technique to the plain diffusion and 
the diffusive reacceleration models as described in 
[14]. We consider here only the effect of the prop- 
agation; the effect of a more realistic source distri- 
bution will be considered elsewhere (Moskalenko 
et al., in preparation). The quality of the fit can 
be judged from Figure 1, which shows the prop- 
agated abundances of elements minus abundances 
measured by ACE. The fitted values are isotopic 
abundances while the plot shows fractional devi- 
ations of elemental abundances for clarity. It can 
be seen that the propagated abundances are repro- 
duced with a maximum error of 5% (10% for the 
least abundant elements). 

The deviations from "0" are mainly due to the er- 
rors in the cross sections. It can be illustrated us- 
ing Carbon as an example. The propagated abun- 
dance of ^^C agrees perfectly with the data since 
the source abundance of this isotope is adjusted it- 
eratively. The isotopic abundance of ^''C cannot be 
adjusted in the same way since it is already zero 
(from the fitting), but it is overproduced due to 
the spallation of heavier nuclei (see discussion in 
[12]). The entire excess of Carbon ^5% is due 
to the overproduction of ^"^C during the CR prop- 
agation. Similarly, F is entirely secondary and the 
excess is due to the cross-section uncertainties. 
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Relative abundance ( Si=100) 
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Figure 2: Source isotopic abundances. Circles: solar system abundances [7], squares: diffusive reaccelera- 
tion model, triangles: plain diffusion model. Isotopes of the same element are connected by the lines. The 
lowest values are included to show the solar abundances; in these cases the CR source abundances are not 
always reliable. 



The isotopic abundances (A > 6) for the two 
models are compared with latest solar abundances 
from [7] in Figure 2. We do not show the iso- 
topes where the source abundances are extremely 
small (<10~^) because the accuracy of the cross 
sections is insufficient in these cases. The derived 
source isotopic abundances generally agree better 
with the latest solar abundances from [7], espe- 
cially the iron group, than with the earlier version 
of solar abundances by Anders and Grevesse [1]. 
Still, the source abundances derived in both models 
are underabundant relative to solar in many cases, 
as is well known, e.g. C, N, O, Ne, S, CI, Ar, Sc. In 
other cases we confirm the remarkable agreement 
with solar: Na, '"'Ca, Mg, Al, Si, ^^Cr, Fe, Co, 
Ni. In many cases a good agreement with solar 
abundances is apparent also for isotopes with large 



proportion of secondary production: Ne, ^Na, 

25,26]y[-g 27^j 29,30(gj 31p 54,57,58pg 59(^q 61jyjj 

The abundances are very much dependent on the 
isotopic production cross sections which are un- 
certain to a large degree, e.g., F, P, Ar, K, Sc, Ti, 
V (see discussion in [12]). The source abundance 
of ^^O, which is mostly primary, is less than so- 
lar by a significant factor. The fragmentation cross 
sections are known much better than the isotopic 
production cross sections, so this result is rather 
robust. Similarly, we confirm that the source abun- 
dances of ^^N and ^°Ne are a factor of 5-6 be- 
low solar The well-known excess ^^Ne/^"Ne (see, 
e.g., [2]) is also evident. Subject to further veri- 
fication and the accuracy of the production cross 
sections, the study shows that many isotopes in CR 
are mostly secondary, ^^C, ^"^O, ^^Ne, 33,34,365^ 

37C1, 38,40^1-, 4lK, 42,43,44ca, 53q. ^his is in ad- 
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dition to Li, Be, B, F, P, Sc, Ti, V that have been 
discussed in the literature for a long time. There 
are also a few cases where the model predictions 
are different from each other, e.g., ^^N, ^^O, ^^Ne, 
^^S, ^^Mn, and some others. These differences can 
be used to constrain or rule out some propagation 
models. 

Interestingly, radioactive ^^Ca and ^-^Mn appear to 
be significantly present in the sources. Radioactive 
^'^Ca is a K-capture isotope that decays (to ^^K) 
with a period T1/2 = 1-03 x 10^ yr and has a negli- 
gible abundance in the solar system. A similar case 
is radioactive ^^Mn, another K-capture isotope that 
decays to ^^Cr with a period T1/2 = 3.74 x 10^ yr. 
Their abundances can be used to constrain the CR 
acceleration time scale. 

Conclusion 

This is the first time that a 'realistic' (i.e. 
fuU spatial- and energy-dependence) propagation 
model has been used to derive isotopic source 
abundances for a full range of CR nuclei. As 
is well-known, the elements with low first ion- 
ization potential (FIP) appear to be more abun- 
dant in CR sources relative to the high-FIP ele- 
ments, when compared with the solar system mate- 
rial. This might imply that the source material in- 
cludes the atmospheres of stars with temperatures 
~10'' K [3]. A strong correlation between FIP and 
volatility (most of low-FIP elements are refractory 
while high-FIP elements are volatile) suggests that 
CRs may also originate in the interstellar dust, pre- 
accelerated by shock waves [10, 5]. CR data tend 
to prefer volatility over FIP, but uncertainties in the 
derived source abundances have prevented an un- 
ambiguous solution. 

The analysis presented in this paper shows that the 
radioactive ''^Ca and ''■'Mn appear to be signifi- 
cantly present in the sources. This may be used 
to study the CR acceleration timescale and seems 
to give more support to the volatility hypothesis. 
The detailed discussion and error analysis will be 
presented in a forthcoming journal pubUcation. 
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